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ABSTRACT 

We use the GALPROP model for cosmic-ray propagation to obtain a new 
estimate of the Galactic component of gamma rays, and show that away from 
the Galactic plane it gives an accurate prediction of the observed EGRET in- 
tensities in the energy range 30 MeV - 50 GeV. On this basis we re-evaluate 
the extragalactic gamma-ray background. We find that for some energies previ- 
ous work underestimated the Galactic contribution at high latitudes and hence 
overestimated the background. Our new background spectrum shows a positive 
curvature similar to that expected for models of the extragalactic emission based 
on the blazar population. 

Subject headings: diffusion — cosmic rays — ISM: general — diffuse radiation 
— gamma rays: observations — gamma rays: theory 
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1. Introduction 

The extragalactic diffuse 7-ray background emission (EGRB) is a superposition of all 
unresolved sources of high energy 7-ray emission in the Universe. Active galactic nuclei 
(AGN) are the dominant class of 7-ray emitters known to emit up to the highest energies. 
There is a consensus that a population of unresolved AGN may contribute to the EGRB, 
however predictions range from 25% up to 100%. Contributions from other extragalactic 
sources have been plausibly suggested: galaxy clusters (Ensslin et al. 1997), energetic parti- 
cles in the shock waves associated with large scale cosmological structure formation (Loeb & 
Waxman 2000; Miniati 2002), or distant 7-ray burst events. Potentially, if reliably derived, 
the EGRB can also provide important information about the phase of baryon-antibaryon 
annihilation (Gao et al. 1990; Dolgov & Silk 1993), evaporation of primordial black holes 
(Hawking 1974; Maki, Mitsui, & Orito 1996), annihilation of so-called weakly interacting 
massive particles (WIMPs) (Jungman, Kamionkowski, & Griest 1996), and extragalactic IR 
and optical photon spectra (Stecker 1999). 

The EGRB is the component of the diffuse emission which is most difficult to determine. 
Its spectrum depends much on the adopted model of the Galactic background which itself is 
not yet firmly established. The isotropic, presumably extragalactic component of the diffuse 
7-ray flux was first discovered by the SAS-2 satellite and confirmed by EGRET (Thompson 
& Fichtel 1982; Sreekumar et al. 1998). However, it is not correct to assume that the isotropic 
component is wholly extragalactic since even at the Galactic poles it is comparable to the 
Galactic contribution from inverse Compton scattering of the Galactic plane photons and 
CMB (Strong, Moskalenko, & Reimer 2000; Moskalenko & Strong 2000). The determination 
of the EGRB is thus model-dependent and influenced by the adopted size of the Galactic 
halo, the electron spectrum there, and the spectrum of low-energy background photons which 
must be derived independently. 

Extensive work has been done (Sreekumar et al. 1998) to derive the spectrum of the 
EGRB based on EGRET data. The relation of modelled Galactic diffuse emission to total 
measured diffuse emission was used to determine the EGRB, as the extrapolation to zero 
Galactic contribution. The derived spectral index —2.10 ± 0.03 appears to be close to that 
of 7-ray blazars. 

In a companion paper (Strong, Moskalenko, & Reimer 2004) we use the GALPROP 
code to infer a new model for Galactic diffuse continuum 7-rays. This model reproduces 
successfully diffuse 7-ray emission from the entire sky and gives a good linear prediction for 
observed vs. predicted 7-ray intensities. In view of the success of this model, we use it as 
the basis of a new determination of the EGRB using EGRET data from 30 MeV to 50 GeV. 
The GALPROP propagation code and previous results are described elsewhere (Moskalenko 
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et al. 2002; Strong et al. 2004, and references therein). 


2. The procedure 
2.1. EGRET data 

The details of the procedure of handling the EGRET data and convolution procedure 
are described in detail elsewhere (Strong et al. 2004); here we provide a brief summary. 

We use the EGRET counts and exposure all-sky maps in Galactic coordinates with 0.5° 
binsize, as in Strong et al. (2000). The sources of the 3EG catalogue have been removed by 
the procedure described in Strong et al. (2000), fully consistent with the 3EG point source 
listings. Apart from the most intense sources, the removal of sources has little influence 
on the comparison with models if sufficiently large sky segments are investigated. For the 
spectra, the statistical errors on the EGRET data points are very small since the regions 
chosen have large solid angle; the systematic error dominates and we have adopted values in 
the range 10-30% depending on energy as shown in Table 2, cf. 15% adopted by Sreekumar 
et al. (1998) and Esposito et al. (1999). The predicted model skymaps are convolved with the 
EGRET point-spread function as described in Strong et al. (2000). Here we use additional 
EGRET data in the energy ranges 10-20, 20-50 and 50-120 GeV. Because the instrumental 
response of EGRET determined at energies above 10 GeV is less certain compared to energies 
below 10 GeV, it is required to account for additional systematic uncertainties. In particular 
the EGRET effective area can only be deduced by extrapolation from the calibrated effective 
area at lower energies (Thompson et al. 1993a). We accordingly adopt values of 0.9, 0.8, and 
0.7 times the 4-10 GeV effective area, respectively. 


2.2. Optimized model for the Galactic diffuse emission 

In the companion paper (Strong et al. 2004) we compared a range of models of Galac- 
tic diffuse emission, based on our CR propagation code GALPROP, with data from the 
Compton Gamma Ray Observatory. There we exploit the fact that the models predict quite 
specific behaviour for different sky regions and this provides a critical test: the “correct” 
model should be consistent with the data in all directions. We show that a new model, 
with moderate changes of electron and nucleon spectra relative to the “conventional” model, 
can well reproduce the 7-ray data and is compatible with locally observed particle spectra 
considering the expected level of spatial fluctuations in the Galaxy. The 7- ray data com- 
parisons were extended over the entire sky and to 100 GeV in energy. We also exploited 
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the recent improved measurements of the local proton, Helium, as well as antiproton, and 
positron spectra which are used as constraints on the proton spectrum in distant regions. 

To fit the Galactic diffuse emission we used 7 test regions covering the sky (for de- 
tails see Strong et al. 2004). The model uses the proton (and He) spectral shape at high 
energies derived from the local data. The secondary antiproton and positron data tracing 
the proton spectrum on a large scale provide an important constraint on the intensity nor- 
malization of the average spectrum. The adopted average electron spectrum resembles the 
local one renormalized upwards by a factor of ~4 and consistent with synchrotron index 
measurements. The adopted electron and nucleon spectra are compatible with the direct 
measurement considering fluctuations due to energy losses and the stochastic character of 
cosmic ray sources and propagation. 

The optimized model fits the observed Galactic diffuse y-ray spectra in all test regions in 
the energy range 30 MeV - 100 GeV. The proposed scenario implies a substantial contribution 
from IC at all energies, but especially below 100 MeV and above 1 GeV. Also IC dominates at 
latitudes |6| > 10° at all energies. The agreement of longitude profiles with the EGRET data 
is generally good considering that the model does not attempt to include details of Galactic 
structure. The agreement of latitude profiles with EGRET is also good, in particular the 
reproduction of the high-latitude variation confirms the importance of the IC component 
which is much broader than the gas-related 7r°-decay and bremsstrahlung emission. The 
outer Galaxy latitude profiles are in excellent agreement with the data. 


2.3. Method 

Since despite its success the model is nevertheless not exact the best approach is to 
fit the observed intensities with a free scaling factor, so that the EGRB is determined as 
the intercept, thus removing any residual uncertainty in the absolute level of the Galactic 
components. This is the same method as in Sreekumar et al. (1998), the difference lies in 
the model. To reduce the effects of Galactic structure, point sources etc. the fits are made 
excluding the plane, using the sky region 360° < l < 0°, 10° < \b\ < 80°; ideally both IC 
and gas-related components would be left free but they are difficult to separate statistically 
at high latitudes, so we make a linear fit to the total IC+7r°-decay-l-bremsstrahlung, with 
the scaling factor and EGRB as parameters. The fit and errors are based on a simple x 2 
analysis, with ( l , b ) bins 360° x 2° to obtain sufficient statistics (at least 10 counts per bin 
were required). For comparison we also made fits to the entire sky; in this case IC and 
gas-related contributions are easily separated, so that fits with both components free can be 
made in addition to fits to the total. The fits in the two regions then give some indication 
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of the model-dependent systematic error in our EGRB estimates. 


3. Determining the EGRB 

The scaling factors determined for the fits in 360° < l < 0°, 10° < |6| < 80°(Table 1) 
reflect the deviations from the model, visible in Fig. 8 of Strong et al. (2004), and are typically 
between 0.85 and 1.15 which is satisfactory considering the EGRET systematic uncertainty 
~15% indicating a good agreement outside the plane. There is only one region (2-4 GeV) 
where the scaling factor is as large as 1.35. For the all-sky fits the 1-2 and 2-4 GeV factors 
are large in the case of the IC related component, probably related to an underestimate 
of the ISRF or of the electron density. In the range 30-50 MeV, the gas-related factor 
is 1.31 which may indicate underestimated bremsstrahlung, but the EGRET data in this 
region are subject to a correction as a result of instrumental calibrations (Thompson et al. 
1993b; see also discussion of the Kniffen factor in Section 3.1 of Strong et al. 2004) and 
thus more uncertain compared to other energies. These are the least reliable ranges of our 
EGRB determination reflecting the discrepancy in the spectrum mentioned in Section 6 of 
Strong et al. (2004). The EGRB is however not very sensitive to the scaling factor, since 
the intercept requires only a small extrapolation to zero Galactic flux. The observed versus 
(fitted) predicted intensities are shown in Fig. 1 (region |5| > 10°) using the parameters in 
Table 1 and Table 2. These plots give an idea of the reliability of the analysis: the relation 
is linear to good accuracy and the EGRB well determined. 

Table 2 presents the EGRB, integrated over each energy bin since this is given directly 
by the fitting. The values from Sreekumar et al. (1998) (their Table 1) were thus integrated 
(assuming an E~ 2 spectrum) to compare with ours. It is seen that the two fitted regions 
(|£>| > 10° and all-sky) give consistent results, indicating that there is no large systematic 
effect; it shows a model-dependent systematic uncertainty of 5-25%. This is comparable to 
the ~15% systematic uncertainty on EGRET effective area. The adopted energy-dependent 
systematic error for EGRET data is also given in the Table. The total error is obtained by 
combining the statistical and EGRET systematic errors in quadrature. 

Fig. 2 shows the spectrum of EGRB as derived in the present work (Table 3). The data 
points of the differential spectrum were obtained from Table 2 (columns 1-3 and column 
4 above 10 GeV) assuming a power-law spectrum within the energy bin. The two E > 10 
GeV points were taken as given in “All-sky, Total” column, while the 20-50 GeV point was 
plotted as a la upper limit. Above 50 GeV the data were insufficient to give a significant 
result by our method. 
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Our estimates lie significantly below those of Sreekumar et al. (1998), in most energy 
ranges. The positive curvature in our EGRB spectrum is interesting and not unexpected for 
a blazar origin (Salamon & Stecker 1998) but in view of the systematic uncertainties should 
not be taken too literally; a similar, less pronounced effect is present in the Sreekumar 
et al. spectrum. Our spectrum is not consistent with a power law. The 0.1 - 10 GeV 
intensity for the region outside the Galactic plane (360° < l < 0°,10° < |£>| < 80°) of 
(11.1 ± 0.1) xl0~ 6 cm -2 sr -1 s -1 can be compared with (14.5 ± 0.5) xlO -6 cm -2 sr -1 s _1 
from Sreekumar et al. (1998). 

Fig. 3 shows the extragalactic X- and 7-ray background, using the compilation by 
Sreekumar et al. (1998) but using our new EGRET values, and also updated COMPTEL 
results (Weidenspointner et al. 2000). 


4. Further checks for systematics due to model 

In order to check the robustness of our EGRB estimates, we compare estimates based on 
data in different hemispheres and four quarter spheres (Table 4) . This tests for the presence 
of apparent anisotropies in the EGRB values, and thus gives an independent estimate of the 
systematic error due to model inadequacies. In order to have sufficient statistics we consider 
only the integral energy range 0.1-10 GeV. The hemisphere deviations are up to 16% from 
the fit to the full \b\ > 10° region. The quarter sphere values deviate by up to 25% from 
the full region. This may reflect essentially different exposures (e.g., southern hemisphere is 
less exposed than northern hemisphere). The deviations are also probably indicating lack 
of symmetry in the observed Galactic 7-ray emission not reflected in the model. However, 
using all regions together, as we have done for our results, will tend to compensate such 
systematics. 


5. Conclusions 

Based on our new optimized model for the diffuse Galactic 7-ray emission, a new EGRB 
spectrum has been derived. It is lower and steeper than found by Sreekumar et al. (1998); 
it is not consistent with a power-law, and shows some positive curvature as expected for an 
origin in blazars. 


We would like to particularly thank David Bertsch for assistance and discussions on 
the subject of the events and instrumental response of the EGRET telescope above 10 GeV 
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Table 1. Scaling factors of model components corresponding to fits in Table 2. 


10° < \b\ < 80°, All-sky 

Energy, MeV total Total IC Gas 


30-50 

1.008 ±0.04 

1.07 

± 

0.017 

0.93 

± 

0.09 

1.31 

± 

0.16 

50-70 

0.864 ±0.02 

0.88 

± 

0.009 

0.66 

± 

0.05 

1.13 

± 

0.06 

70-100 

0.872 ± 0.01 

0.84 

± 

0.005 

0.82 

± 

0.03 

0.86 

± 

0.02 

100-150 

0.851 ± 0.01 

0.82 

± 

0.004 

0.85 

± 

0.03 

0.80 

± 

0.01 

150-300 

0.874 ± 0.01 

0.83 

± 

0.003 

0.95 

± 

0.02 

0.79 


0.01 

300-500 

0.965 ± 0.01 

0.90 

± 

0.004 

1.14 

di 

0.04 

0.84 

± 

0.01 

500-1000 

1.02 ±0.01 

1.09 

± 

0.004 

1.28 

± 

0.04 

0.87 

± 

0.01 

1000-2000 

1.15 ±0.02 

1.22 

± 

0.007 

1.50 

± 

0.07 

0.98 

± 

0.02 

2000-4000 

1.35 ± 0.04 

1.39 

± 

0.014 

1.54 

dfc 

0.11 

1-11 

± 

0.04 

4000-10000 

0.96 ± 0.07 

0.92 

± 

0.022 

0.98 

± 

0.15 

0.89 

± 

0.07 

10000-20000 

0.98 ± 0.04 

0.81 

± 

0.058 

0.62 

± 

0.39 

0.91 

± 

0.21 

20000-50000 

... 

0.87 

’± 

0.16 

1.14 

:± 

1.3 

0.71 

± 

0.73 


Note. — “Total” refers to fitting sum of all model diffuse components, 
“IC” and “gas” refer to fits separating IC from gas-related components. Only 
statistical errors are given. 




Table 2. Estimates of EGRB obtained by fitting optimized model 500190 to EGRET 

data. 


Energy, MeV 

10° < |&| < 80°, 
total 

Adopted relative 
systematic error 

All-sky 

Total IC±gas 

Sreekumar et al. (1998) a 

30-50 

16.8 ±0.66 

0.30 

16.0 ±0.48 

17.0 ±0.84 

24.0 ± 7.0 

50-70 

10.6 ±0.19 

0.15 

10.4 ±0.13 

11.1 ±0.22 

13.26 ±2.58 

70-100 

6.66 ±0.10 

0.10 

6.7 ±0.065 

6.75 ±0.11 

7.83 ± 1.05 

100-150 

4.48 ±0.07 

0.10 

4.6 ±0.045 

4.55 ±0.076 

5.5 ±0.75 

150-300 

3.92 ± 0.06 

0.10 

4.2 ± 0.040 

3.92 ±0.067 

5.4 ±0.72 

300-500 

1.20 ±0.04 

0.10 

1.39 ±0.025 

1.17 ±0.041 

1.97 ±0.268 

500-1000 

0.76 ±0.04 

0.10 

0.93 ±0.023 

0.70 ± 0.037 

1.36 ±0.185 

1000-2000 

0.32 ± 0.03 

0.10 

0.39 ±0.018 

0.26 ±0.028 

0.617 ±0.084 

2000-4000 

0.24 ± 0.02 

0.10 

0.30 ±0.014 

0.25 ±0.022 

0.304 ± 0.044 

4000-10000 

0.167 ±0.02 

0.10 

0.13 ±0.010 

0.12 ±0.016 

0.1956 ±0.0288 

10000-20000 

0.025 ±0.04 

0.25 

0.034 ±0.011 

0.04 ±0.02 

0.053 ±0.016 

20000-50000 


0.25 

0.011 ±0.022 


0.0276 ±0.0096 

50000-120000 





0.0147 ±0.0063 


Note. — “Total” refers to fitting sum of all model diffuse components, “IC” and “gas” refer to fits separating 
IC from gas-related components. For our fits statistical errors are given together with adopted systematic errors 
on the EGRET response. Units: 10“ 6 cm“ 2 sr _1 s~ x . 

a Values are from their Table 1 integrated over the energy bin. The error bars include systematic errors. 
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Table 3. EGRB intensity. 


Intensity, 

J3 7 , MeV cm" 2 sr*" 1 s” 1 MeV -1 

Total error 

30-50 

8.40 

x 10~ 7 

2.54 

X 

10' 

-7 

50-70 

5.30 

x 10~ 7 

0.80 

X 

10' 

-7 

70-100 

2.22 

x 10 -7 

0.22 

X 

10' 

-7 

100-150 

8.96 

X 

o 

00 

0.91 

X 

10' 

-8 

150-300 

2.61 

x 10~ 8 

0.26 

X 

10' 

-8 

300-500 

6.00 

X 10~ 9 

0.63 

X 

10' 

-9 

500-1000 

1.52 

x 10~ 9 

0.17 

X 

10' 

-9 

1000-2000 

3.20 

x 10~ 10 

0.44 

X 

io- 

-10 

2000-4000 

1.20 

x 10~ 10 

0.16 

X 

io- 

-10 

4000-10000 

2.78 

x 10~ n 

0.43 

X 

io- 

-11 

10000-20000 

3.40 

x 10~ 12 

1.39 

X 

io- 

•12 

20000-50000 

< 1.11 x 10~ 12 


•• 

• 


Note. — EGRB intensity (as plotted in Fig. 2) derived 
from columns 1-3 of Table 2 and column 4 above 10 GeV. 


Table 4. Test for model-dependent systematics. 

1 

b 

Intensity 0.1-10 GeV 

Description 

0-360 

< -10, > +10 

11.10 + 0.12 

N+S hemispheres 

0-360 

< -10 

11.70 + 0.15 

N hemisphere 

0-360 

>+10 

9.28 ± 0.21 

S hemisphere 

270-90 

< -10, > +10 

11.90 + 0.17 

Inner Galaxy N-f S 

90-270 

< -10, > +10 

9.75 ± 0.17 

Outer Galaxy N+S 

0-180 

< -10, > +10 

10.80 + 0.17 

Positive longitudes N+S 

180-360 

< -10, > +10 

11.60 ± 0.16 

Negative longitude N+S 

270-90 

> +10 

13.00 + 0.22 

Inner Galaxy N 

270-90 

< -10 

9.14 + 0.32 

Inner Galaxy S 

90-270 

> +10 

10.60 ± 0.22 

Outer Galaxy N 

90-270 

< -10 

8.18 + 0.34 

Outer Galaxy S 


Note. — Units: 10 6 cm -2 sr -1 s 1 . Errors are statistical only. 






Fig. 1. — Observed versus predicted intensities for region 360° < l < 0°,10° < \b\ < 80 
(GALPROP model ID 500190). Predicted is Galactic only, intercept is EGRB. 
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Ey, MeV 

Fig. 2. — Comparison of our EGRB spectrum (solid, red) as given in Table 3 with that from 
Sreekumar et al. (1998) (dots, magenta). The dashed (green) line is Sreekumar et al. fit: 
2.743 x 10 -3 £ -21 cm -2 s -1 sr -1 MeV -1 . 



Photon Energy (keV) 


Fig. 3. — Extragalactic X-ray and 7 -ray spectrum. Data compilation from Sreekumar et al. 
(1998) except for COMPTEL (Weidenspointner et al. 2000) and EGRET 30 MeV - 20 GeV 
(this work). 





